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^ ; Abstract. 

, The Main Injector Neutrino Oscillation search (MINOS) uses two detectors separated by 735 km 

■ to measure a beam of neutrinos created by the Neutrinos at the Main Injector (NuMI) facility at 
^ I Fermi National Accelerator Laboratory. The experiment has recently reported an observation | ll] 

I . of disappearance consistent with neutrino oscillations. We describe the manner in which the 

^ ' experiment's results depend on the correct understanding and modeling of hadronic systems. 
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^ : INTRODUCTION 

o ; 

X I The Main Injector Neutrino Oscillation Search (MINOS) is a long baseline, two-detector 

^ ■ neutrino oscillation experiment that will use a muon neutrino beam produced by the 

O-i! Neutrinos at the Main Injector (NuMI) facility at Fermi National Accelerator Laboratory 

^ I (FNAL) [0, [l]. The measurement is conducted by two functionally identical detectors, 

• located at two sites, the Near Detector (ND) at FNAL and the Far Detector (FD) in the 

• r-j ! Soudan Underground Laboratory in Minnesota. The experiment began collecting beam 

rN I data in March 2005 and has recently reported disappearance consistent with quasi- 

_c3 \ two-neutrino oscillations according to jArngjI = 2.74+q 26 ^ 

10-3ev2/c4 and sin 26 > 

0.87at68%CL[l]. 

The MINOS detectors are tracking-sampling calorimeters, optimised to measure neu- 
trino interactions in the energy range 1 < ^ 50GeV. The active medium comprises 
4.1 cm- wide, 1.0 cm-thick plastic scintillator strips arranged side by side into planes. 
Each scintillator plane is encased within aluminum sheets to form a light-tight module 
and then mounted on a steel absorber plate. The detectors are composed of a series of 
these steel-scintillator planes hung vertically at a 5.94cm pitch with successive planes 
rotated by 90° to measure the three dimensional event topology. Wavelength-shifting 
and clear optical fibers transport scintillation light from each strip to Hamamatsu multi- 
anode photomultiplier tubes that reside in light-tight boxes alongside the detector. Both 
detectors are magnetized so as to measure muon charge- sign and momentum via curva- 
ture. 

The disappearance measurement is done by using the event topology to identify 
interactions as v^— CC, rather than NC/Ve— CC, then reconstructing the neutrino energy 
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FIGURE 1. events measured in the MINOS calibration detector. The pions are incident from the 
left and two (orthogonal) views of the event are shown. Cells demarcate the detector's scintillator strips, 
with shading/color indicating the pulse-height of hits. The event on the left was collected at 1 Ge V/c beam 
momentum while the one on the right is from a 2 GeV/c run. To remove hits due to optical cross-talk in the 
multi-anode PMTs, only strips registering a pulse-height larger than 1 .5 PE (photoelectrons) are shown. 



as the sum of the muon energy and the energy transferred to the nucleus : Ey = E^ + v 
where nu = y x Ey The oscillation hypothesis is tested by comparing the measured neu- 
trino energy spectrum to the spectrum expected in the absence of oscillations. The latter 
spectrum is anchored to observations made with the Near Detector. MINOS endeavors 
to measure Am^ and sin^20 with an accuracy of better than 10%. The results depend 
upon a reliable knowledge of the event selection efficiency and the energy scale, both 
of which are affected by the modeling of the hadronization process and the detector's 
signal and topological response to those hadrons. We will discuss the way that MINOS 
has dealt with these issues. In addition, uncertainties in the production of hadrons in the 
NuMI target result in uncertainties in the neutrino flux and poor agreement with data. 
We close by describing the manner in which data from the Near Detector was used to 
constrain the simulation of the NuMI beam. 
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FIGURE 2. Pion and electron line shapes as measured in the MINOS calibration detector 10, IH]. The 
Monte Carlo calculation is that of the experiment's GEANT3/GCAL0R simulation. The same program is 
used to model neutrino induced interactions in the MINOS Near and Far Detectors. 



HADRONIC ENERGY SCALE 

MINOS data suggests that (assuming quasi-two-neutrino oscillations) the largest sup- 
pression of occurs for I < Ey < 2GeV and that the suppression slowly shrinks as 
the energy increases to about 10 GeV. This is a relatively low energy range, but not low 
enough to make the experiment dominated by quasi-elastic — CC for which Ey may 
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FIGURE 3. MINOS detector response to K"^ induced showers iQ]. Data measured in the MINOS cali- 
bration detector are compared with the result of GEANT3/GCAL0R and GEANT3GHEISHA simulations. 
The TT+Z response asymmetry is absent above 4 GeV. 



be reconstructed from only the muon energy and scattering angle. Instead, events with 
1 < < lOGeV have, after experimental cuts, an average inelasticity of ^ 0.3. We are 
therefore interested in hadronic showers in the few GeV range and in which the energy 
is carried by a relatively few particles. The response of calorimeters to these low energy 
showers is notoriously difficult to model. Experiments generally remedy this by mea- 
suring the response to single-particles, either to perform a direct calibration or as a way 
of constraining/tuning the shower model. Since the MINOS Near and Far Detectors are 
large and have been constructed underground, direct exposure to a calibration beam was 
not possible. Instead a dedicated Calibration Detector (CalDet) was built to establish the 
energy scale and develop the calibration technique. 

The CalDet, a scaled-down but functionally equivalent model of the MINOS Far and 
Near detectors, was exposed to test beams in the CERN PS East Area during 2001-3 to 
establish the response of the MINOS calorimeters to hadrons, electrons and muons in the 
momentum range 0.2-10 GeV/c. The detector consisted of sixty 2.5 cm x 1 m x 1 m steel 
plates each of which supported a plane of twenty-four 100 cm-long scintillator strips. 
Alternating planes were rotated ±90° as in the larger MINOS detectors. The sampling 
fraction is 6.4% and the detector has a longitudinal (transverse) granularity of O.lSAj 
(0.25 Ai). Identification of ;r + e and p was accomplished using a time of flight system 
and several threshold Cerenkov counters. Pions and muons were discriminated using the 
event topology. Fig. [H shows two pion induced showers measured in the detector. Data 
were collected in 200MeV/c steps (1 GeV/c for Pbeam > 4GeV/c) in both positive and 
negative beam polarity and in two beamlines (Tl 1 and T7). 

The detector was calibrated using a procedure similar to that employed in the Near 
and Far detectors [@]. Gain and light-output non-uniformities were corrected using an 
LED based light injection system [7] and through-going cosmic-ray muons [8]. The 
absolute scale, in "muon equivalent units", was defined as the average signal per- 



plane induced by muons with momenta between 0.5-1.1 GeV/c at the time of plane 
crossing [9]. The momentum was determined by requiring that the muons stop in the 
detector (p^ < 2.2GeV/c), then working backward along the track, inverting the range- 
energy relationlO. 

Figure [21 shows and electron line-shapes measured in the calibration detector 
and compared to the GEANT3/GCAL0R [ill, HES] Monte Carlo. The simulated re- 
sponse to electrons was in good agreement with the data after upstream energy loss 
(e.g., via bremsstrahlung and ionization) was included in the simulation. The e/7t re- 
sponse ratio is energy dependent but always > 1 and ^1.27 for momenta above a 
few GeV. The line-shapes display an asymmetric, high-side tail characteristic of non- 
compensating calorimeters. The agreement between the data and the MC is surprisingly 
good, even at the lowest momentum settings. The energy resolution may be parame- 
terised as 56%/ Ve © 2% for 7t^ showers and 21 .4%/v^ © 4% /E for electrons. 

The detector's response to is shown as a function of energy in Fig. [31 The data 
are in somewhat better agreement with the GCALOR based simulation than they are with 
the GHEISHA [Il4l [l5ll based one. The agreement is excellent for induced showers 
but poorer for showers. This reflects a measured asymmetry, as much as 6% for few- 
GeV energies, in the detector's response to and n^. This asymettry was not present 
in the response to nor in the range of /i^ and vanishes for Ej,; >3 GeV. Its origin 
may lie in differences, between -|-A^ and +N, in the multiplicity of secondaries 
produced in inelastic interactions at these low energies [16]. 

Hadronic showers produced in neutrino interactions are simulated with the version of 
GEANT/GCALOR shown here. Because these showers consist of multiple particles, and 
because the response to e, ;r and p all differ (the latter two below ~ 1 GeV where protons 
tend to range out before interacting hadronically) we cannot derive the relation between 
detector signal and energy carried by the shower particles directly from the CalDet data. 
Instead we derive the relationship from the MC but assign a somewhat pessimistic 6% 
systematic error on the response. This error is small enough to have a negligible effect on 
the precision of results from the low-statistics first year MINOS dataset. The uncertainty 
may be reduced by introducing data based correction factors into the relation or, more 
ideally, by improving the hadronic shower model. Over the long term MINOS plans to 
transition to a GEANT4 based simulation in the hopes of improving the hadronic shower 
treatment. 

FINAL STATE INTERACTIONS 

Test beam measurements, such as the one discussed above, are used to determine a de- 
tector's hadronic and electromagnetic response to single particle induced showers. The 
single particle based calibration may be used for energy measurements of multi-particle 
showers (e.g., "jets") if the energy response is sufficiently linear and if the particle con- 
tent of the shower is well known or, especially, if the calorimeter's response to different 
particles of the same energy is identical (e.g., a compensating calorimeter) [17]. These 
concerns apply to neutrino induced showers (which bear some similarity to jets) but 
there is an additional complication owing to the fact that the neutrino interacts in the 
dense nuclear medium of the target nucleus. Hadrons produced by the neutrino must 
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FIGURE 4. Left: Energy distribution of simulated V^-CC events at the MINOS Far Detector. The 
effect of neutrino oscillations is included for Am^ = 2.47, 2.72 and 3.12 x lO^^eV^ (68% CL limits of 
Ref. Ol]) and sin^20 = 1. Individual curves correspond to uncorrected ±1(7 variations in the formation 
time, pion absorption cross-section and the neutrino induced 2 ;r production cross-section. The MC sample 
corresponds to the "high-statistics limit", approximately x 10 of the planned MINOS exposure. Right: 
The variation in the calorimeter's response as a function of the energy transfer V and under different 
final state interaction scenarios. Zero on the vertical scale corresponds to the n npnp model of pion 
absorption and To=0.52fm/c. 



escape the nucleus to create a signal in the detector but, at MINOS energies, often suffer 
a hadronic interaction before doing so. It is important to account for these "final state 
interactions" because neutrino oscillation experiments seek to measure the energy, v, 
transferred by the neutrino to the target, so as to reconstruct Ey =E^ + v. The presence 
of final state interactions causes v 7^ E^i^, where the latter is the visible energy of shower 
particles exiting the nucleus. We will describe the way in which final state interaction 
are modelled and how uncertainties are accounted for. 

Neutrino interactions in MINOS are modelled by NEUGEN-v3 [Hi]. NEUGEN in- 
cludes a final state interactions (also referred to as "intranuclear re-scattering") package 
known as INTRANUKE. The code is anchored to a comparison of final states mv -\-d 
and V + Ne interactions as measured in the BEBC and ANL-12 ft bubble chambers fl^. 
The library includes a treatment of pion elastic and inelastic scattering, single charge 
exchange and absorption in a cascade simulation of the final state. The relative proba- 
bilities for these processes are approximately 35:50:10:5 (40:30:7:23) at a pion kinetic 
energy of 1 GeV (250 MeV). The formation zone concept is included by suppressing 
interactions that would have occurred before the pion has travelled a distance I = Top/m 
To = 0.52fm/c [20]. The simulation includes a treatment of pion absorption inspired 
by [12 in . Absorbed pions transfer their energy to a npnp cluster and the individual nucle- 
ons are then tracked in GEANT. A more comprehensive description of the INTRANUKE 
code and it's effect in MINOS is given in Ref. [22]. 

We derive the effect that parametric uncertainties in the final state interaction model 
have on the oscillation result as follows. First, the model contains a number of param- 
eters for which the uncertainty may be estimated, however roughly, based on external 
data. The dependence of the reconstructed shower energy on v was studied as a function 



of these model parameters and the most influential were identified. The parameters were 
then varied in an uncorrelated way over multiple trials and the visible — CC energy 
spectra were constructed, accounting for the detector energy smearing. This was done 
without oscillations so as to represent measurements made in the Near Detector and with 
oscillations at varying Am^ and sin^20 to represent measurements at the Far Detector. 
An example of simulated Far Detector spectra is shown on the left in Fig. ID The Am^ 
values were chosen according to the best-fit and ±lc7 uncertainties on the initial MINOS 
result lHI]. For each trial a fit was performed to extract Am^ and sin^20 and the differ- 
ences with respect to the input values were histogrammed. The mean and width of the 
distributions is an estimate of the uncertainty attributed to the parameter uncertainties in 
the final state interaction model. We find 5Am^ ^3x 10^^ eV^/c^ and 5 sin^ 20 ^ 0.01. 
These errors are small compared to the overall uncertainty as may be noticed from the 
fact that the curves for different values of Am^ in Fig. |4] are easily distinguished. 

A much larger uncertainty occurs as a consequence of the rather simple way in 
which the final state interaction code simulates pion absorption and inelastic collisions. 
Absorption is modelled, by default, as solely due to ;r — npnp but we know that 
processes such as ;r ^ npp and n nnp may also occur and that the nucleons are 
themselves attenuated inside the nuclear medium (this effect is not accounted for). 
In addition, n + N inelastic collisions are modelled but the energy lost by the n is 
discarded and secondary n are not created. These model based deficiencies result in 
a large uncertainty on the amount of visible energy emitted from the nucleus. 

The model based uncertainties cannot be accounted for by the parametric variation 
technique described above. Instead, we study simulations done (a) without re-scattering, 
(b) with pion absorption according to ;r — npnp (the default) and (c) in which the energy 
of absorbed pions is discarded. The differences a — b and c — b in the reconstructed 
shower energy are shown as a function of v on the right-hand side of Fig. IH The curve 
labelled "OFF" corresponds to a — c and represents the change in visible energy due to 
inelastic pion interactions. The curve labelled "Full n absorption" corresponds to c — b 
and represents the change in visible energy due to pion absorption. For comparison the 
(much smaller) effect of ±lc7 variations in the formation time is also shown. To estimate 
the uncertainty on the oscillation parameters, the v ^ £^shower relation is modified 
according to Fig. |4] and the analysis is repeated. We find that the parameters shift by 
dAnP' 6 X 10~^ eV^/c^ and 5 sin^ 20 ~ 0.05. When drawing the oscillation contours, 
uncertainties in both the single particle energy scale and final state interactions are 
included as an 1 1% penalty term ("nuisance parameter") on the hadronic energy scale. 

The experiment is currently working on improvements in the final state interaction 
model [23] to remedy the deficiencies mentioned above. In addition, the model is 
being re-implemented to facilitate re-weighting of MC events. A re-weightable model is 
extraordinarily important because it vastly reduces the number of MC events needed 
to estimate systematic uncertainties. Re- weighting machinery would also allow the 
parametric variation technique described above to be used directly in the oscillation fit 
by including the final state interaction model parameters as nuisances. The importance 
of event by event re-weighting is one reason that experiments may prefer homegrown 
codes, for which the model parameters can be altered and uncertainties understood, over 
more sophisticated but also more opaque general purpose codes. 
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FIGURE 5. The reconstructed energy distribution of V^— CC events observed in MINOS. The data 
were collected in the (a) low (b) medium and (c) high energy beam configurations and are compared to 
the predictions of the neutrino beam Monte Carlo, using FLUKAOS to simulate hadron production in the 
NuMI target. We show the MC prediction before and after the tuning procedure described in the text. 



ESTIMATION OF THE NEUTRINO FLUX 

The neutrino beamline model is separated into three parts: (a) a simulation of the hadrons 
produced by 120GeV/c protons incident on the NuMI target and (b) the propagation 
of those hadrons and their progeny through the magnetic focusing elements, along the 
677 m decay pipe, and into the primary beam absorber and (c) decay of the mesons and 
calculation of the probability that any neutrino progeny traverses the Near and Far detec- 
tors. Hadronic production is simulated using a detailed model of the target geometry and 
material composition. By default, interactions are generated with FLUKAO 5 [ 24ll25l . l2^ 



but output from the MARS-vl5 [|27i S |29J and GEANT-FLUKA [30] codes, as well 
as calculations based on the parameterisations of BMPT 113 ill and Malensek [|32ll . are 
used as cross-checks. The produced hadrons are propagated in a GEANT3 simulation of 
the NuMI beamline. Decays in which a neutrino is produced are saved and later used as 
input for neutrino event simulation in the Near and Far Detectors. 

The neutrino fluxes predicted by the different models vary by 10-20%. This is un- 
suprising given the relatively paucity of hadron production data in the Xf,pT region 
interesting to NuMI, not to mention at 120GeV/c and on a thick Carbon target f33|]. 
Therefore, the rather good agreement between the v^ — CC energy spectrum measured 
in the Near Detector and the MC prediction based on FLUKAOS (Fig. [S]), comes as a 
pleasant suprise. The data/MC agreement is, however, not perfect and can be improved 
by taking advantage of the flexibility of the NuMI beamline. The NuMI target may be re- 
motely moved along the beam axis, changing the Xf , pr region focused by the magnetic 
horns, and thereby the beam energy. Figure [5] shows data taken in three target postions: 



z — 10cm, 100 cm and 250 cm where z=Ocm is the locatio of the target when fully in- 
serted into the first horn. These data, along with data taken at z = 10cm with the horns 
off and with the horn current modified by ±8%, are used in a parameteric fit to constrain 
the beam model. 

The fit is based on a parametric model fsT], similar to that of BMPT, which is able to 
describe the d^N /dxfdpT distribution predicted by FLUKA. The model parameters were 
allowed to float in a fit to the v^— CC energy spectra. The fit also included nuisance pa- 
rameters to account for uncertainties in beam focusing, intensity, and detector response. 
The mean transverse momentum of pions, (pj) was constrained to the FLUKA0 5 value 
of 364MeV/c, with a penalty term of 15 MeV/c, obtained from the variation of hadron 
production models. The results of the fit are shown as the "Full MC Tuning" curve in 
Fig. [51 The fitting procedure dramatically improves the data/MC agreement, particularly 
in the lowest energy configuration in which most of the oscillation dataset was recorded. 
The oscillation analyses use these results as the first step in the prediction of the neutrino 
spectrum at the Far Detector. 

SUMMARY 

We have attempted to highlight some of the ways in which MINOS depends upon the 
correct understanding and modelling of hadronic interactions from 120 GeV (hadron 
production) to a few tens of MeV and less (calorimetry). Hadronic interactions are diffi- 
cult to model correctly and whenever possible we have attempted to benchmark and/or 
constrain simulations using internal or external datasets. One expects that models will 
never completely agree with the data so much effort is concerned with quantifying the 
effect that disagreements have upon the oscillation results. Uncertainties are generally 
included as nuisance parameters and fitting is greatly facilitated by the ability to re- 
weight events for model changes. Finally, in the interest of brevity we have been forced 
to omit some important topics, such as the effect of shower modeling uncertainties on the 
NC background estimated in the reconstructed v^— CC sample as well as the sensitivity 
of the Vg— CC event selection to the fragmentation/hadronisation treatment. 

The author would like to thank J. Morfin, H. Gallagher, R. Ransome and S. Dytman 
for advice and comments regarding final state interactions and N. Mokhov for comments 
about low energy pion nucleus interactions. The hadron production fitting is due to 
S. Kopp, Z. Pavlovic and R Vahle. 
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